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ABSTRACT 

The 2:1 mean-motion resonance with Jupiter harbours two distinct groups of 
asteroids. The short-lived population is known to be a transient group sustained in 
steady state by the Yarkovsky semimajor axis drift. The long-lived asteroids, however, 
can exhibit dynamical lifetimes comparable to 4 Gyr. They reside near two isolated 
islands of the phase space denoted A and B, with an uneven population ratio B/Aczt 
10. The orbits of A-island asteroids are predominantly highly inclined, compared to 
island B. The size-frequency distribution is steep but the orbital distribution lacks 
any evidence of a collisional cluster. These observational constraints are somewhat 
puzzling and therefore the origin of the long-lived asteroids has not been explained so 
far. 

With the aim to provide a viable explanation, we first update the resonant popula¬ 
tion and revisit its physical properties. Using an IV-body model with seven planets and 
the Yarkovsky effect included, we demonstrate that the dynamical depletion of island 
A is faster, in comparison with island B. Then we investigate (i) the survivability of 
primordial resonant asteroids and (ii) capture of the population during planetary mi¬ 
gration, following a recently described scenario with an escaping fifth giant planet and 
a jumping-Jupiter instability. We also model the collisional evolution of the resonant 
population over past 4 Gyr. Our conclusion is that the long-lived group was created 
by resonant capture from a narrow part of hypothetical outer main-belt family during 
planetary migration. Primordial asteroids surviving the migration were probably not 
numerous enough to substantially contribute to the observed population. 

Key words: minor planets, asteroids: general - methods: numerical. 


1 INTRODUCTION 

The 2:1 mean-motion resonance with Jupiter, hereinafter 
denoted as J2/1, is one of the major first-order Jovian res¬ 
onances intersecting the main asteroid belt. It is associated 
with the Kirkwood ga p known as the Hecuba gap llKirkwoodl 
1 18671 : fSchwcizcr 1969|), which is commonly considered to be 
the borderline separating the outer main belt and the Cy- 
bele region. Its name is derived from (108) Hecuba which is 
an outer main-belt asteroid but technically does not exhibit 
the exact 2:1 commensurability (ISchubartll 19641) . 

The 2:1 resonance, which was originally believed to 
be depleted of asteroids, harbours several small-sized bod¬ 
ies as it was realised by observations in the 20th century 
(the first resonant asteroid (1362) Griqua was recognized 
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in 1943 H abel[l9 59)). The follow-up studie s based on an¬ 
alytical and semi-analytical method s (e.g. iMurravl Il986l : 
iMoo ns. Morbkielh & Migljorinil 11998T) . symplectic mapping 
fe.g.lRoig fe Ferraz-Mellc jl999l) and frequency analysis (e.g. 
Nesvorny & Ferraz -Mcllo 19971 1 helped to get insight to the 
internal structure of the resonance and to investigate the 
character of various resonant orbits affected by a compli- 
cate d interplay between the J 2/1 and overlapping secu- 
l ar (iMorbidelli fe Moons [l993|)_and secon dar y re sonances 
l|Wisdomll 19871 : Henrard. Watanabe fe Moon J 19951 ') . In par¬ 
ticular, it was discovered that while the regions of the 
overlap give rise to a strongly chaotic behaviour of the 
asteroidal orbits, there are two separate regions of quasi¬ 
regular motion l ocat ed inside t he reso n ant part of the phase 
space (|Franklinl 1 19941 : iMichtchenko Ferraz-Mellol Il997l ). 
These regions were designated as stable islands A and B 
dNesvornv fe Ferraz-Mellall997l . see e.g. Fig. [l]). 
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However ,_it_ was _also_suggested by 

iFerraz-Mello. Michtchenko fe Roiel (Il998l l that the sta¬ 
bility of the islands could have been weakened in the past. 
The mechanism responsible for this increase of chaos was 
introduced as 1:1 resonance of the J2/1 libration period 
with the period of Jupiter-Saturn great inequality (GI). 
Although the present value of the GI period is 880 yr, 
it was probably smaller w hen configuration of plan etary 
orbits was more compact (IMorbidelli fe Cridal l2007fl and 
it slowly increased during Jupiter’s and Saturn’s divergent 
migration. If the value of the GI period was temporarily 
comparable with the libration period of the J2/1 asteroids 
(which is typically ~ 420 yr), the aforementioned resonance 
would have occurred. 

With ongoing improvements of the observation tech¬ 
niques and sky surveys, there is now a population of aster¬ 
oids inside the J2/1 which is large enough to be analyzed sta¬ 
tistically with numerical methods. Long-term integrations in 
a si mplified model with fo ur gian t planets were performed by 
iRoig. Nesvornv fc Ferraz-Mellol (2002) who identified both 
short-lived and long-lived asteroids within the observed pop¬ 
ulation, which consisted of 53 bodies by then. Investigat¬ 
ing resident lifetimes of the resonant asteroids, IRoig et ~ahl 
(l2002ll demonstrated that part of the population escapes 
from the resonance on time-scales ~10Myr, while other as¬ 
teroids may have dynamical lifetimes comparable to the age 
of the Solar System. The orbits of these long-lived aster¬ 
oids were embodied in island B and part of them was also 
present in its boundaries, being affected by chaotic diffusion. 
Surprisingly, the stable island A was empty. 

In a series of papers fo cuse d on the 2:1 mean-motion res - 
onance, iBroz etahl ( 2005 ) and iBroz fc Vokrouhlickil (j2008l f 
identified significantly more bodies as resonant. The latter 
catalogue contained 92 short-lived and 182 long-lived as¬ 
teroids. They also realised that nine of the long-lived as¬ 
teroids were located in island A, while the rest was resid¬ 
ing in the island B and its vicinity. They successfully ex¬ 
plained the origin of the short-lived asteroids by a numer¬ 
ical steady-state model, in which the resonant population 
was replenished by an inflow of outer main-belt asteroids 
driven by the Yarkovsky semimajor axis drift. However, this 
model fails in the case of the long-lived population which was 
thought to be created by Yarkovsky induced injection of the 
nearby Themis family asteroids that exhibit simi lar incli¬ 
nations as the B-island asteroids. iBroz et aD (|2005l J showed 
that the objects transported from Themis are usually per¬ 
turbed shortly after entering the resonance and rarely reach 
the islands. 

[Roig et ahl (|2002l l argued that the steep size-frequency 
distribution (SFD) of the long- lived asteroids, which was 
inconsistent with a steady state dDohnanvilll969li . could im¬ 
ply recent collisional origin. On the contrary, there was no 
collisional cluster identified in the orbital distribution by 
iBroz fe Vokrouhlickvl (120081 '). thus a collisional family would 
have to be older than 1 Gyr. Moreover, the presence of the 
highly inclined A-island asteroids would require unrealistic 
ejection velocities, assuming that the disruption occurred in 
the more populated island B. 

As there is no self-consistent explanation, the main goal 
of this paper is to provide a reasonable model for the forma¬ 
tion and evolution of the long-lived population, taking into 


account the paucity of bodies in island A, the differences in 
the inclination and the non-equilibrium size distribution. 

The long dynamical lifetimes of the dynamically sta¬ 
ble asteroids, together with the failure of the aforemen¬ 
tioned hypotheses, strongly indicates that their origin may 
be traced back to the epoch of planetary migration. Re¬ 
cent advances in migration theories suggest that a primor¬ 
dial compact configuration of planetary orbits and subse¬ 
quent violent planetary migration might have led to a desta¬ 
bilization of several regions that are otherwise stable un¬ 
der the current plane tary configuration. T he pop ulatio ns o f 
Trojans and Hildas iN esvornv . Vokrouhlickv fe Morbidellil 


120131 ; IRoig fe Nesvornvll2014 ) may serve as examples of such 
a process: the primordial populations in the 1:1 and 3:2 reso¬ 
nances with Jupiter were totally dispersed and the observed 
populations were formed by resonant capture of bodies, orig¬ 
inating either in the outer main belt or in the transneptunian 
disc of comets. It seems inevitable that the 2:1 resonance 
with Jupiter also undergoes significant changes of its loca¬ 
tion, inner secular structure and asteroid population, during 
the planetary migration. 

The structure of the paper is as follows. We first use the 
latest observational data from the databases AstOrb, Ast- 
DyS, WISE and SDSS to update the observed population 
in Section [5] We also briefly review the characterization of 
resonant orbits and describe a method of dynamical map¬ 
ping. In Section [3] we study whether the planetary migra¬ 
tion may cause depletion or repopulation of the long-lived 
J2/1 asteroids. The dynamical models we create are based 
on simulations with prescribed evolution of planets in con¬ 
text of the modern migration scena rio with five giant p l an¬ 
ets a nd jumping-Jupiter instability (Ne svornv & Morbid clli 
|2012|J . Section []] is focused on dynamical simulations cover¬ 
ing the late stage of planetary migration during which the 
GI period evolution might have influenced the stability of 
the long-lived asteroids. In Section[S] we examine the effects 
of collisions on the long-lived population. Finally, Section [6] 
is devoted to conclusions. 


2 OBSERVED RESONANT POPULATION 

In this section, we first describe methods we use for identi¬ 
fication and description of resonant orbits as well as for dy¬ 
namical mapping of the resonance. O ur approach is similar 
to the one used bv iRoig et alJ d2002l l. We identify resonant 
orbits in a recent catalogue of osculating orbital elements 
and we study their dynamical lifetimes on the basis of long¬ 
term numerical integrations in a simplified model with four 
giant planets only. Our goal is also to survey available data 
for physical properties of the population, namely the abso¬ 
lute magnitudes, sizes and albedos. In the last subsection, 
we revisit results of Skoulidou et al. (in preparation) in or¬ 
der to compare the simplified four-giant-planet model with 
a more sophisticated framework including terrestrial planets 
and the Yarkovsky effect. Especially, we look for the differ¬ 
ences in the dynamical decay rates in the stable islands A 
and B. 
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2.1 Characterization of resonant orbits 

The 2:1 resonance critical angle is defined as: 

a — 2Aj — A — vj , (1) 

where Aj and A are the mean longitudes of Jupiter and of 
an asteroid, respectively, vj is the asteroid’s longitude of 
perihelion. The critical argument of any body trapped inside 
the J2/1 resonance librates (quasi-periodically changes on 
an interval < 27 t) with a typical period of about 420 yr. On 
the other hand, the critical argument of the asteroids outside 
the resonance circulates. This provides us a useful tool for 
the identification of the resonant asteroids. 

The libration of a is linked to periodic changes of the 
osculating semimajor axis a, the eccentricity e and the incli¬ 
nation I. These changes are coupled together as described 
by the adiabatic invariant N of the asteroid’s motion in the 
circular and planar restricted Sun Jupiter-asteroid system 

N = \Ja ^2 — y/l — e 2 cos . (2) 

The presence of other planets and variable eccentricity of 
Jupiter’s orbit give rise to multiple perturbations in the J2/1 
and prevent integrability of the orbits. 

None the less, the a, e, I coupling is preserved to a cer¬ 
tain degree. Semimajor axis oscillates around the libration 
centre, which is positioned approximately at the exact res¬ 
onance a res ~ 3.27 AU. The eccentricity and inclination at¬ 
tain their maximal values when the oscillation of a is at its 
minimum and vice versa. 

An inconvenient consequence of this behaviour is that 
the standard av eraging methods for th e computation of 
proper elements dKnezevic fc ]VIilanill2003h do not retain any 
information about the libration amplitude (i.e. the proper 
semimajor axes of all resonant asteroids approach the value 
of ares)- Therefore the proper elements are not the appro¬ 
priate choice when studying resonant orbits. 

In order to properly characterize the libration ampli¬ 
tude, we use an alternative set o f resonant (or pseudo¬ 
proper) elements llRoig et al.l 12002 .1. The idea is to record 
the osculating orbital elements at the moment when they 
reach their extremal values during the libration cycle. These 
values can be found as the intersections with a suitably de¬ 
fined reference plane in the osculating elements space. A set 
of conditions determining this plane can be written as 

a = 0 A > 0 A vj — vj j = 0 A fl — flj = 0 , (3) 

where f i denotes the longitude of node and subscript J is 
used for Jupiter. The purpose of the conditions for vj and fl 
is to eliminate secular variations of the resonant elements. 
The whole set enables the resonant elements to be recorded 
when the osculating semimajor axis a reaches its minimum, 
the eccentricity e and the inclination I attain maxima. Note 
that as a consequence, the J2/1 asteroids are always depicted 
on the left-hand side (closer to the Sun) of the libration 
centre in the resonant elements space (see e.g. Fig. 0 . 

Moreover, temporal evolution of resonant elements may 
serve as the first indicator of the stability. The reason is that 
stable orbits exhibit stable libration with only a small vari¬ 
ation around the mean value. Therefore, successive intersec¬ 
tions with the reference plane do not move significantly and 
the recorded resonant elements are nearly exact constants 


of motion. The situation for the unstable orbits is just the 
opposite and the intersections slowly disperse in time. This 
fact propagates numerically to the resulting resonant ele¬ 
ments, if we compute their standard deviation as an error 
of the mean value obtained over a significant period of time 
(~100 kyr). 

Nevertheless, due to higher order perturbations and sec¬ 
ular effects, the above conditions are seldom satisfied exactly 
and one has to use less confined criteria in the following form 
when numerically integrating the orbits: 

|oj < 5° A > 0 A \vj — - CC7J | < 5° . (4) 

We use the difference between successive numerical time 
steps (denoted as A a) rather then the time derivative. Note 
that we completely omit the condition for Q. This simplifi¬ 
cation can be compensated by verifying that the maximal 
value of the inclination is reached when recording the reso¬ 
nant elements. 


2.2 Dynamical maps 


In this section, we summarize our approach to dynamical 
mapping of the 2:1 resonance and its secular structure. Our 
method is verified by comparing the map computed for the 
present configuration of planets with th e separatrices and lo - 
cations of secular resonances derived bv lMoons et al.l lll998li . 
An important outcome is our ability to record global changes 
inside the resonance when the planetary orbits are recon¬ 
figured and also to examine dynamical stability in various 
regions of the resonant phase space. 

Our method employs the definition of the resonant ele¬ 
ments. The respective values should change systematically in 
the case of unstable resonant orbits, whereas stable orbits 
should only oscillate with small variations. The same be¬ 
haviour is exhibited by actions of a dynamical system; time 
series of their extremal or mean v alues are the r efore often 
used for d ynamical mapping (e.g Laskar j 19941 : iMorbideiil 
1 19971 : iTsiganis. Knezevid fe Varvoglisll2007l h 

We proceed as follows. First, we divide an investigated 
part of the phase space into the boxes of the same size. Let 
the central coordinates (a, e, I) of each box represent a set of 
initial resonant elements (this is accomplished by setting the 
osculating angular elements so that condition (0 is fulfilled). 
We integrate the initial orbits over several millions of years 
and track the evolution of the resonant elements. Finally, 
we determine the differences <5a r , 5e r and <5sin/ r between 
the initial and final (or the last recorded) resonant elements. 
These values represent the dynamical stability of the initial 
orbit and we use them as a characteristics of the entire box. 

For the purpose of expressing the total displacement in 
the phase space, we define the distance 


d = 



(<5e r ) 2 + (5 sin J r ) 2 , 


(5) 


where a r denotes the arithmetic mean of the initial and 
final resonant semimajor axis. The distance d is similar 
to the metric used in the Hierarchical Clustering Method 
dZappala et al.lll995|j for family identification. 

In order to verify our method, we constructed a dynam¬ 
ical map of the 2:1 resonance in the present configuration 
of planets. We investigated the phase space in the intervals 
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Figure 1. A dynamical map of the 2:1 mean-motion resonance 
with Jupiter computed in the present configuration of giant plan¬ 
ets. Top panel: the map is plotted in the resonant semimajor axis 
a r vs resonant eccentricity e r plane and it is averaged over five sec¬ 
tions in the resonant inclination sin I r . The color coding of boxes 
represents the average distance in the phase space travelled by a 
test particle with initial orbital elements placed within the box. 
Additionally, separatrices and borders of the secular resonances, 
which were computed bv lMoons et al.|(|l998h . are plotted over the 
map. The left solid line corresponds to the separatrix of the J2/1. 
The near-vertical solid line is the libration centre of the J2/1. The 
dotted lines indicate the and z /5 resonance, and the dashed 
line represents the Kozai resonance. We mark the stable islands 
denoted A and B and also the region of overlapping secondary 
resonances (Roig & Ferraz-Mcllo 1999). Bottom panel: the map 
plotted in the resonant semimajor axis a T vs sine of resonant in¬ 
clination sin 7 r plane. It was computed for test particles with the 
resonant eccentricity e r = 0.25. 

a G (3.195,3.275) AU, e G (0.1, 0.5), I G (0°,25°) and split 
them into a grid of 40 x 40 x 5 boxes. With the aim to im¬ 
prove our statistics, we randomly generated two more test 
particles in the vicinity of each central particle (the disper¬ 
sion of particles is not exceeding 20 per cent of the box size). 
I 11 this way we created the initial conditions for 24, 000 test 
particles. 

In the next step, we integrated the orbits for £ S pan = 
10 Myr. For all t hese simulations, we use the symplectic 
integrator swift (jLevison fe Duncanl 1 19941) with a built - 
in second-order symplectic scheme (jLaskar Sz Robutel|[200lh 
and with an implementation of digital filters for the com¬ 


putation of the resonant elements based on criterion © 
(Br oz et al .l l2005f ). This symplectic integrator enables us to 
use a time-step At = 91.3125 d. The simulations we perform 
are simplified: as we study the outer main-belt we take into 
account only the gravitational interactions with the Sun and 
the four giant planets. The terrestrial planets are neglected, 
except for a barycentric correction applied to the initial con¬ 
ditions. We use the Laplace plane to define our initial frame 
of reference. We do not consider any non-gravitational ac¬ 
celeration at this point. 

The final step was to calculate the distance d for each 
test particle. The average value d for the particles initially 
placed in the same box was taken as a measure of the dy¬ 
namical instability in the box. Note that there is no temporal 
measure in this method, i.e. the map reflects the maximal 
displacement in each box during the integration timespan 
tspan but it does not say how fast did this displacement oc¬ 
cur. 

We plot the mean values d in the top panel of Figure |T| 
The map is an average projection of all sections in sin J r to 
the (a r ,e r ) plane. It can be co mp ared with several separa¬ 
trices derived by [Moons et al] (il998i 'l. It is clear that the 
map reflects the inner secular structure of the 2:1 resonance 
very well: we can locate both the stable islands A and B, 
separated by the secular resonance. The borderline at 
higher values of e r is formed by the Kozai and v 5 separa¬ 
trices. The low-eccentricity region near the libration cen- 
tre is affected by presence of multiple secondary resonances 
(|Roig fc Ferraz-Mellolll999 ). 

The bottom panel of Fig. [T] displays the projection of 
our map to the (a r ,sin/ r ) plane but this time for a small 
interval of e r because the shape of the separatrices and stable 
islands strongly depends on e r . The map demonstrates how 
the shape of the stable islands changes with the resonant 
inclination J r . 

Let us conclude that the suitability of the resonant ele¬ 
ments for the dynamical mapping serves as an independent 
confirmation that they indeed reflect the regularity of res¬ 
onant orbits and they retain important properties of the 
proper elements at the same time. 

2.3 Lifetimes of asteroids in the 2:1 resonance 

Let us now discuss our approach to the identification and 
classification of the resonant asteroids. We numerically prop¬ 
agated the orbits of known numbered and multi-opposition 
asteroids in the broad surroundings of the J2/1 to identify 
those trapped inside, using a time series of the critical argu¬ 
ment a. We extracted the osculating ele ments o f the main- 
belt objects from the AstOrb database dBowe U l2012ll as of 
2012 November 15 to set-up the initial conditions for the first 
short-term integration. Our choice of the borders in the os¬ 
culating ( a,e ) plane was ei = 0.45 (a — 3.24) / (3.1 — 3.24) 
and e2 = 0.5 (a — 3.24) / (3.46 — 3.24). We ended up with 
11,469 orbits. The corre spon ding planetary epliemeris were 
taken from JPL DE405 dStandishll2004l f for the given Julian 
date. We numerically integrated the orbits for 10 kyr and 
we recorded the critical argument a for each asteroid during 
the simulation. We have found 374 librating asteroids. 

The following long-term integration requires a different 
approach. Our goal is to determine the future orbital evolu¬ 
tion of the resonant asteroids and estimate their dynamical 
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lifetime. However, we are limited by strong chaotic diffusion 
in the J2/1 resonance. In other words, even a slight change 
of the initial conditions can significantly alter the orbital 
evolution. Because the orbital elements of each observed as¬ 
teroid are only known with a finite accuracy, one should 
consider all values within the observational uncertainty as 
possible initial conditions to cover all alternatives of future 
orbital evolution. 

To account for this chaotic behaviour, we first matched 
the AstOrb data for librat ing asteroids with correspond¬ 
ing AstDyS uncertainties dKnezevic fe Milanil 120031 ). The 
matching failed in four case^U only and hence we discarded 
these bodies from the J2/1 population. We then used a 
pseudo-random generator to create a bundle of 10 synthetic 
orbits for each asteroid which are close to its nominal orbit. 
The generated orbital elements fall within the Gaussian dis¬ 
tribution (over ±3 ct interval) in the nonsingular osculating 
element space. These synthetic test particles are called close 
clones. In this manner, we obtained 4,070 orbits (1 nom¬ 
inal and 10 synthetic for each body) which we integrated 
over 1 Gyr. By this procedure, we can study several possible 
realisations of future orb ital motion. 

Following again iRoig et ~akl J2002I) . we define the dy¬ 
namical lifetime r as the asteroid’s timespan of residence 
inside the resonance. Test particles leaving the resonance 
are usually discarded due to highly eccentric or inclined or¬ 
bits, which lead to planetary crossing or fall into the Sun. We 
calculated the median dynamical lifetime t as the median 
value of the residence lifetimes of close clones. 

We divided the J2/1 asteroids into three groups by 
virtue of their median dynamical lifetime f as follows: 

• f ^ 70 Myr: short-lived/unstable/Zulus. Number of 
identified bodies: 140. 

• f G (70,1000) Myr: long-lived/marginally sta- 
ble/Griquas. Number of identified bodies: 106. 

• f ^ 1 Gyr: long-lived/stable/Zhongguos. Number of 
identified bodies: 124. 

For reference, the numb ers of resonant asteroids identified 
in the previous paper fBroz fe Vokrouhlickvl d2008l ) were 92 
short-lived and 182 long-lived asteroids. 

We used the same long-term integration of the resonant 
orbits to construct orbital distribution of asteroids in the 
J2/1. We calculated average values of the resonant elements 
for each asteroid and each close clone over the time interval 
of 1 Myr; then we computed the arithmetic mean for the 
asteroid and its close clones together. We define the standard 
deviation of the resonant elements as the uncertainty of the 
arithmetic meaifl- 

We depict the resulting mean resonant elements as pro- 

1 These particular bodies have the following designations: 
2003 E031, 2003 HH4, 2006 SE403 and 2006 UF152. The rea¬ 
son for the databases mismatch might be in different methods 
used for orbit identification. 

2 Let us note that the standard deviations of the mean resonant 
elements will be slightly overestimated. The reason for this arises 
from our approach to the generation of initial conditions because 
the initial elements of the close clones are generated randomly 
thus resulting in uncorrelated sets of quantities even though they 
should be correlated, as described by the correlation matrices 
available in orbital elements databases. 



resonant semimajor axis a r [AU] 



Figure 2. Orbital distribution of the J2/1 asteroids in the reso¬ 
nant semimajor axis a T vs the resonant eccentricity e r plane (top) 
and in a r vs sine of the resonant inclination sin I T plane (bottom). 
The symbols correspond to the dynamical lifetime of each body: 
filled circles denote dynamically stable Zhongguos, empty circles 
marginally stable Griquas, crosses unstable Zulus. The error bars 
indicate standard deviations of the computed orbits. We show the 
same borders and structures as in Fig. [l] Only a relatively small 
part of the unstable population is depicted. 


jections to the (a r ,e r ) and (a r ,sin/ r ) planes in Fig. [5] The 
orbital distribution exhibits similar properties as in the pre¬ 
vious studies. Zhongguos reside at the very centre of the 
stable islands, 11 Zhongguos in island A and 113 in island 
B. The majority of B-island Zhongguos have low resonant in¬ 
clinations while A-island Zhongguos are more inclined. Gri¬ 
quas partially overlap with the Zhongguo group, but they 
either reside on orbits with lower semimajor axis or higher 
inclination than Zhongguos. There is no clear separation be¬ 
tween the orbits of Griquas and Zhongguos. It is therefore 
questionable whether the division of the long-lived popula¬ 
tion into these two groups has a solid physical foundation, 
apart from the dynamical lifetime. The orbital position of 
Griquas with respect to Zhongguos indicates that Griquas 
might be chaotically diffusing part of Zhongguos which have 
haphazardly entered the periphery regions of the stable is¬ 
lands. 

Since we focus on the long-lived population, a great part 
of the short-lived bodies is not displayed in Fig. [2] (they are 
located outside the range of the plot). The depicted short- 
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geometric albedo p v 


Figure 3. A histogram of the albedo distribution among the J2/1 
asteroids: the visual geometric albedo py vs the n umber of aster¬ 
oids N . The data for 44 asteroids were available in lMasiero et al] 

1 I 201 ih . 


lived asteroids are chaotically drifting towards the Kozai 
separatrix. 


2.4 Albedo, colour and size-frequency 
distributions 


In order to study the basic physical properties of the J2/1 
population, we first searched the WISE database. We were 
able to extract the visual geometric albedos pv and the effec- 
tive diameters D , inferred from NEATM thermal models by 
iMasiero et al.l (l201ll ). for 44 asteroids inside the J2/1 (out 
of 370). 

The data obtained allow us to examine the albedo dis¬ 
tribution in the resonance (Fig. [3j. The lowest albedo is 
(0.034 ± 0.002) while the largest value is (0.24 ± 0.03). This 
relatively large value corresponds to asteroid 2001 RN2 and 
indicates that this object likely belongs to the S taxonomic 
type. The majority of the asteroids (98 per cent) have albedo 
lower or equal to (0.136 ± 0.004). The shape of the albedo 
distribution is typical for the outer main-belt region where 
C-types dominate fe.g. IPeMeo fe Carrvll2013l '). 

As the next step, we investigated the diameters of as¬ 
teroids. For the 44 cases we simply use the value D and its 
uncertainty. For the rest of the asteroids we calculat e their 
app roximate diameter using the following relation (jHarrisI 
1199811 : 


D,^10- 


( 6 ) 


where we insert the absolute magnitude H from the AstOrb 
catalogue and we use the mean albedo pv = (0.08 ± 0.03). 
The standard deviation of D is then computed as the prop¬ 
agated uncertainty of both H and pv- 

We construct the cumulative size-frequency distribu¬ 
tions (SFDs) of the resonant population and individual 
groups, as shown in Fig. [4] We fit the steep part of the 
distribution with a power-law function N (> D) oc D 1 to 
estimate the slope parameter 7 for further comparison. The 
value of 7 clearly depends on the chosen interval of diam¬ 
eters over which we approximate the SFDs with the power 
law. We set the nominal fitted interval as D £ (7.5,18) km. 
The SFDs start to bend at the lower limit of this nominal 


group 

nominal 7 

variation of 7 

short-lived 

- 3.2 

(- 2 . 5 ,- 3 . 7 ) 

long-lived 

- 4.3 

(-3.7,-5.1) 

Griquas 

- 3.0 

(-3.0,-3.3) 

Zhongguos 

- 5.1 

(-3.9, -5.1) 


Table 1 . The slopes 7 resulting from the method of least squares 
which we applied to fit the SFDs of dynamical groups with a 
power law function .Y (> D) oc DY The middle column shows 
the result of the fit in the nominal interval of diameters D 6 
(7.5,18) km, the third column reflects the variation of 7 if the 
limits of the interval are shifted. 



diameter D [km] 


Figure 4. Size-frequency distributions (SFDs) of individual dy¬ 
namical groups residing inside the J2/1: the diameter D vs the 
cumulative number N (> D) of asteroids larger than D. A loga¬ 
rithmic scaling is used in the plot. The distributions correspond 
to the short-lived population (top left) and long-lived popula¬ 
tion (top right), which is then divided to marginally stable Gri- 
quas (bottom left) and stable Zhongguos (bottom right). The 
steep part of each distribution is approximated with a power-law 
N (> D ) oc DY The nominal borders of the fitted regions (rang¬ 
ing from 7.5 km to 18 km) are shown by vertical lines. The value 
of the slope 7 is also given. We also plot stationary Dohnanyi-like 
slope (7 = —2.5) for comparison. 


interval except for SFD of Griquas which remains steep up 
to D ~ 4 km. We also record the variation of the slope 7 by 
slightly changing the fitted range of diameters. The results 
are su mmarized in Table [T] Unlike the results of lBroz et al.l 
d2005l l. the SFD of updated Griquas is not s hallo wer than 
a Dohnanyi-like SFD with 7 = —2.5 llDohnanvil 1 19691 ) . In 
Section 0 we shall use the observed SFDs for comparison 
with the outcomes of our collisional models. 

Finally, we als o searched the SDSS MOC catalogue 
(lParker et al.ll2008l ). We found astrometric and photomet¬ 
ric data for 81 resonant bodies for which we constructed 
a colour-colour diagram (see Fig. \B- The principal compo¬ 
nent denoted as a* is defined on the basis of measurements 
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Figure 5. The optimized colour a* vs the i — z colour diagram 
of the J2/1 resonant asteroids found in the SDSS MOC cata¬ 
logue. Error bars represent the standard deviations of the dis¬ 
played quantities. The outliers with a* > 0 are labeled with cor¬ 
responding catalogue numbers. 

in filters r, i and g as a* = 0.89 (g — r) + 0.45(r — i) — 0.57 
and it can be used to distinguish C-complex (a* < 0) and 
S-complex (a* > 0) asteroids. The distribution in the dia¬ 
gram is typical for outer main-belt C-type asteroids, but one 
can also see several outliers with a* > 0. We did not found 
any significant relation between the colours and the orbital 
distribution. 


2.5 The long-term stability of the islands A and B 

The observed resonant population has likely evolved 
over Gyr-long time-scales. Hence, it is important 
to assess its dynamical stability over similarly long 
time-spans. This pr o blem w a s rec ently revisited by 
ISkoulidou. Tsiganis fc Varvoglisl (120141 '). Here we will use 
the main results of this study that are relevant for our 
work; a complete dynamical study will be presented in a 
different paper (Skoulidou et al. in preparation). 

When studying the outer asteroid belt, we typically 
ignore the weak perturbative effects of the terrestrial 
planets, as thi s spee ds up the numerical propagation. 
ISkoulidou et al.l J2014ri showed that this is not a safe op¬ 
tion, when studying the 2:1 resonant population for time- 
spans longer than 1 Gyr, as numerous weak encounters with 
Mars actually induce small-scale chaotic variations. These 
variations build-up slowly with time and eventually assist in 
breaking the phase-protection mechanism that would other¬ 
wise prevent asteroids on moderately eccentric orbits (e ~ 
0.4) from encountering Jupiter on this time-span. Moreover, 
when the Yarkovsky effect is added, asteroids with D < 
20 km can escape from the re sonance more efficient l y. Thi s 
is because, as also shown in iBroz fe Vokrouhlickvl (|2008h . 
the combined action of resonance and slow Yarkovsky drift 
forces the asteroid orbits to slowly develop higher eccentric¬ 
ities as a consequen ce of adiabatic invarian ce. The compu¬ 
tations presented in ISkoulidou et all d2014l'l show that, in a 
physical model that takes into account the aforementioned 
phenomena, the 2:1 population would decay roughly expo¬ 
nentially in time, with an e-folding time of order 1 Gyr. 



Figure 6. Time evolution of the number of resonant objects, 
(TV (t) /No), normalized to the initial value. This is the result of 
a 3 Gyr-long simulation, within the framework of the 7-planets 
model and including Yarkovsky-induced drift in a. The two 
curves refer to the two islands (A=open triangles and B=solid 
circles). Superimposed are the least-square fitted exponentials, 
yielding e-folding times of t\ = (0.57 ± 0.02) Gyr and rg = 
(0.94 ± 0.02) Gyr, respectively. 

In Fig. [S] we present the results of a similar simula¬ 
tion that spanned 3 Gyr, in a model that contained the 
seven major planets (Venus to Neptune; ‘7pl’ model). This 
simu lation treate d th e 2:1 population, as was defined by 
IBroz fc Vokrouhlickvl ll2008l l. Also, an approximate treat¬ 
ment of the Yarkovsky acceleration was included in the equa¬ 
tions of motion, assuming a constant drift rate in semimajor 
axis equal to 2.7 x 10 -4 D -1 AUMyW 1 . For each object its 
diameter D was estimated, using its catalogued H value and 
assuming an albedo of 0.06-0.08, and a value of either 0° or 
180° obliquity was randomly assigned. The population was 
found to decay exponentially in time, according to what was 
described in the previous paragraph. 

However, what is more interesting here is that the two 
sub-populations that are contained in the two quasi-stable 
resonant islands are not diffusing out at the same rate. Fig. [6] 
shows two sets of curves (one for each island): clearly island 
A is depleted faster than island B. Fitting exponentials on 
the data, we get the corresponding e-folding times: ta = 
(0.57 ± 0.02) Gyr for island A and r B = (0.94 ± 0.02) Gyr 
for island B. This uneven depletion is an important dynami¬ 
cal property of the resonance and is likely one of the reasons 
behind the observed A/B asymmetry. The importance of 
this observation will become more evident in the following 
sections. 


© 2015 RAS, MNRAS POQ.fTlESOl 








































































8 O. Chrenko et al. 


3 EFFECTS OF THE JUMPING-JUPITER 
INSTABILITY 

We investigate the role of the jumping-Jupiter instability on 
the depletion and eventual repopulation of the J2/1 in this 
section. We study survival of a hypothetical primordial pop¬ 
ulation, and resonant capture from t he outer main belt, both 
in the fifth giant planet scenario jNesvornv fc Morbidellil 
l2012ll . Our simulations cover only the instability phase itself, 
demonstrating at first which process is plausible. 


3.1 Simulations with prescribed migration 

We adopt the technique suggested bv iNesvornv et all d2013l l 
for simulating the jumping-Ju piter instability. We mod ified 
the SWIFT_rmvs3 integrator (|Levison fc Duncanlll994l '). so 
that the evolution of massive bodies is given by a prescribed 
input file and the evolution of test particles is computed by 
the standard symplectic algorithm. This method allows us 
to exclude the transneptunian disc of planetesimals from our 
integrations, because its damping and scattering effects are 
no longer needed; planets evolve ‘the way they are told’ by 
the input file. Moreover, the exact progress of the migration 
scenario is then exactly reproducible. 

The evolution of migrating planets is prescribed 
by the fif th giant plane t scenario, developed in 
iNesvornv fe Morbidellil (I2012T ) (see Fig. [7]). This scenario 
has already been proved reliable in terms of reproducing the 
orbital architecture of planetary orbits, the period ratios 
and secular frequencies, and several distributions of minor 
solar-system bodies. The prescribed input file is the result 
of a complete self-consistent simulation with five giant 
planets and a massive transneptunian disc of planetesimals 
(A/disc = 20M Ea rth)- We use only a 10 Myr portion of 
the simulation, containing the jumping-Jupiter phase. The 
input data sampling is Ati npu t = 1 yr which is sufficient for 
our purpose. 

Our integrator transforms the input data to the Carte¬ 
sian coordinates (i.e. orbital elements to positions and ve¬ 
locities) and interpolates them according to the integration 
time step At = 0.25 yr. The interpolation is done by a for¬ 
ward/backward drift along Keplerian ellipses, starting at 
input data point preceding/forthcoming to required time 
value. The results of both drifts are averaged using the 
weighted mean, where the weight depends on temporal dis¬ 
tance between the data points and the required time value. 


3.2 Capture from the outer main belt 

We investigated a possibility that the resonant asteroids 
were captured from outer main-belt orbits near the J2/1 
during Jupiter’s jump. We present a simulation of the cap¬ 
ture in this section. 

Initial conditions. Initially, we distributed 5, 000 test par¬ 
ticles uniformly over the region which is about to be covered 
by the moving 2:1 resonance. This can be estimated easily by 
Kepler’s third law since we know the evolution of Jupiter’s 
semimajor axis a priori. Our choice of the semimajor axes 
distribution was a £ (3.06, 3.32) AU. The distribution of ec¬ 
centricities and inclinations was also chosen as uniform and 
corresponds to the extent of a moderately excited main belt: 



Figure 7. Orbital evolution of giant planets i n the f i fth gi- 
ant planet scenario, adopted from lNesvornv &; Morbidellil 820121 1. 
during the jumping-Jupiter instability, as it was reproduced by 
our modified integrator. We plot the time t vs the semimajor axis 
a, the pericentre q and the apocentre Q. Each evolutionary track 
is labeled with the name of the corresponding giant planet. 


e £ (0,0.35), I £ (0°, 15°). The angular elements were ran¬ 
domly distributed over (0°,360°) interval. 

Evolution of test particles. We recorded the time series 
of resonant elements of the test particles during the integra¬ 
tion and process ed t h em off -line using the Savitzky-Golay 
smoothing filter (IPress et al 120071) with 0.1 Myr range of the 
running window and a second order smoothing polynomial. 

The result of our simulation is shown in Fig. [H] Note 
that for the sake of simplicity, we use the resonant elements 
for all particles, even for those not trapped inside the J2/fl 
Shortly after the onset of migration, the test particles still 
retain the uniform character of the initial distribution. A 
few minor mean-motion resonances can be identified. The 
apparent gap between the resonance centre and the synthetic 
population is intentional, because we do not want any of the 
test particles to be initially placed inside the resonance. 

After 2 Myr, the relaxation processes start to take place. 
The 2:1 resonance changes its position with respect to the 
initial one and starts to perturb several eccentric orbits at 
the outskirts of the synthetic population. All of the per¬ 
turbed as teroids fall i nto t he short-lived unstable zone (dis¬ 
cussed in iBroz et~ah ][ 2 QQl). The minor mean-motion reso¬ 
nances begin to pump the eccentricity of the test particles 
residing inside them. 

At 6.26 Myr, an instantaneous discontinuity, known as 
Jupiter’s jump , occurs. The 2:1 resonance, co-moving with 
jumping Jupiter, changes its location suddenly and a large 
number of test particles flow through the libration zone, 
while others are excited and ejected, or left behind in the 
emerging inner Cybele region. 

Finally at 10 Myr, a population between the left-hand 
separatrix and the libration centre is stabilised, as the result 


3 The resonant elements of the non-resonant asteroids do not 
have any special physical meaning, but they correspond to the 
extremal secular variations of osculating elements as given by 
equation 0. 
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resonant semimajor axis a r [All] 


Figure 8. A simulation of the resonant capture from the outer main belt in the fifth giant planet scenario iNesvornv &; ATo r liiclel 1 il 2(J 1 2i) . 
We plot the evolution of test particles in the resonant semimajor axis a T vs resonant eccentricity e r plane. The vertical line indicates an 
approximate location of the 2:1 resonance centre which migrates inwards together with Jupiter. The integration time is shown in three 
of panels (t = 0, 2, 6.26 and lOMyr) and corresponds to the timeline in Fig. [T] The dashed arrows and lines indicate an approximate 
extent of the libration zone in case of the migrating and stable resonance, respectively. From top to bottom and left to right, the plots 
show: the initial conditions, the relaxed population of test particles prior to jump, the state during the jumping-Jupiter instability and 
the final state. Note that all test particles are depicted in terms of resonant elements for simplicity (even the non-resonant orbits). 


of the resonant capture. The bodies outside the left-hand 
separatrix are simply remnants of the initial population of 
test particles. A few bodies between the libration centre and 
the right-hand separatrix are in fact outside the 2:1 reso¬ 
nance, as well as the asteroids in the Cybele region. 

Captured population. Instead of another long-term inte¬ 
gration, we use the rather efficient dynamical mapping intro¬ 
duced in Section EU to check the stability of the captured 
population in the post-migration configuration of planets. 
This allows us to identify the counterparts of the stable is¬ 
lands A and B in the phase space for this planetary configu¬ 
ration. Hence, we can trace the long-lived candidates of our 
captured bodies simply by selecting particles with resonant 
elements falling within the range of the islands. 

The dynamical map for the 2:1 resonance is shown in 
the Fig. [21 It covers the phase space in the intervals a G 
(3.115,3.195) AU, e G (0.1,0.5) and I G (0°,25°), which 
were divided into a grid of 40 x 40 x 5 boxes (and we finally 
took an average over all sections in the inclination). We set 
up and followed three test particles per box for up to 10 Myr. 

Comparing the map with Fig.Q] one can see very similar 
structures. The stable islands A and B are there, separated 
by the secular resonance. The Kozai resonance separa¬ 
trix can also be easily identified. A major difference is the 
shape and size of the A island - the resonance is not 


present in the depicted part of the phase space, which ef¬ 
fectively enlarges the A island with respect to the present 
state. The selected candidates for long-lived orbits are also 
shown in Fig. [9] Note that there are very few A-island bodies 
captured on highly eccentric orbits, therefore the aforemen¬ 
tioned difference in shape of island A should not strongly 
affect the results. Using the dynamical mapping, we identi¬ 
fied N^ nth = 69 candidates in island A and N^ nth = 254 
candidates in island B. 

The next step is a rescaling of the initial population so 
its particle density would reach realistic values. We used 
the observed particle density in the main belt for bod¬ 
ies with diameters D ^ 5 km within the following region: 
a G (2.95, 3.21) AU, e G (0, 0.35), / G (0°, 15°). The intervals 
of e and I correspond to those of our synthetic initial pop¬ 
ulation. On the other hand, the interval of semimajor axes 
was shifted into the region which is currently not depleted 
by the presence of the J2/1. We further increased this par - 
ticle density by a factor of three (IMinton fe Malhotrall20ld) 
to account for the dynamical depletion of the main belt after 
the reconfiguration of planets during the last ~ 3.85 Gyr. If 
we increase the particle density throughout the initial popu¬ 
lation, the number of captured long-lived bodies larger than 
5 km would be N s ^ led = 1552 and A^ caled = 5857 in island 
A and B, respectively. 

Finally, we assume that the dynamical depletion rate 
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ta = (0.57 ± 0.05) Gyr and tb = (0.94 ± 0.05) Gyr. We 
use t — (3.9 ±0.1) Gyr as a value of the decay time pe¬ 
riod. Combining the results together and identifying the 
limit values, our model predicts that we should observe 
^ymodei _ l on g_li ve ci bodies larger than 5km in island 
A and jV^ 101161 = 62-121 asteroids in island B. The observed 
values are N X bs = 2 and N B bs = 71. The results of our model 
and the observed values are in good agreement, which sup¬ 
ports the resonant capture scenario. 

Turning our attention to the bottom panel in Fig. [9] we 
can discuss whether the inclination distribution of captured 
synthetic bodies can evolve towards the observed one. It can 
be clearly seen that the number of low-7 asteroids in island 
A is lower than the number of high-7 asteroids. Assuming 
strong depletion of the whole island and considering that 
only 1-3 objects larger than 5 km may survive up to the 
present according to our model, the low-7 asteroids are more 
likely to be depleted completely. That is in agreement with 
the observations because the observed asteroids in island A 
reside exclusively on highly inclined orbits. 

On the other hand, the population captured in island B 
stretches almost uniformly over a relatively large interval of 
inclinations. Considering the shape of the real island shown 
in Fig.Q] which shrinks with increasing inclination, the long¬ 
term diffusion in the high-7 region will partially deplete this 
part of the captured population. But in order to reproduce 
the observed state in island B, it is necessary for the source 
population to contain low-7 asteroids predominantly. 


3.3 Survival of the primordial population 


Figure 9. A final result of resonant capture from the outer main 
belt in the fifth giant planet scenario. Top: the resonant semi¬ 
major axis a r vs the resonant eccentricity e r . The background of 
the plot is the dynamical map of the 2:1 mean-motion resonance 
with Jupiter computed for the post-migration configuration of the 
giant planets. All white symbols indicate orbits of test particles 
captured during the simulation (they correspond to the final state 
in Fig. [Sj but only a part of the phase space containing the stable 
islands is plotted). Test particles embedded in the dark regions 
of the dynamical map are candidates for a long-term stability 
and we distinguish them by circles. Other particles are marked 
by crosses. Bottom: the resonant semimajor axis a T vs sine of the 
resonant inclination sin 7 r . We plot only the particles captured 
inside the stable islands. Black circles indicate test particles cap¬ 
tured in island B and gray open squares indicate those captured 
in A island. The horizontal dashed line is plotted for reference, 
because vast majority of the observed B-island Zhongguos reside 
on low inclinations sin7 r ^0.1 and all of the observed A-island 
Zhongguos have high inclinations sin7 r 0.1. 

due to chaotic diffusion is the same for the post-migration 
islands as for the islands in the observed configuration. 
Knowing the approximate number of captured long-lived as¬ 
teroids, we employ an exponential decay law and the ex¬ 
tremal valuesQ of the e-folding times in island A and B 


Let us study hypothetical long-lived primordial orbits and 
their survival during the planetary migration. The proce¬ 
dure is similar to that in the previous section, the major 
difference is of course in the initial conditions setup. Here 
we aim to study only the bodies which reside inside the J2/1 
when the instability simulation begins. Another important 
property that should be satisfied is the long-term stability 
of the initial orbits - otherwise one could unintentionally 
simulate different effects such as survival of the short-lived 
population, etc. 

Initial conditions. How to create a set of test particles 
on the long-lived orbits at the beginning of the migration 
scenario? The only viable solution is to use the dynamical 
mapping once again. To this effect, we slightly modify the 
standard procedure described in Section m First, the in¬ 
tegration is not done in a stable unvarying configuration of 
planets. Is is instead carried out during the first 5 Myr of the 
prescribed migration scenario, starting in the pre-migration 
configuration of planets and evolving into a configuration 
closely preceding the instability. 

Second, we use only the recorded changes of the res¬ 
onant inclination <5sin7 r when constructing the map (i.e. 
we omit Sa T and Se T when evaluating the metric given by 
equation (f5jl ). This is necessary because as Jupiter migrates 


4 These values were derived in Section 12.51 but we use a slightly 
higher standard deviation (0.05 instead of 0.02). The reason is 
that the e-folding times were obtained from a model with seven 
planets and the Yarkovsky drift included, with long-lived aster¬ 


oids of all sizes. In this section, however, we focus on asteroids 
larger than 5 km and consequently, the e-folding times are rather 
upper limits for us. We thus artificially increase the original de¬ 
viation to account for this difference. 
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resonant semimajor axis a r [AU] 

Figure 10. A dynamical map of the 2:1 mean-motion resonance 
with Jupiter computed during the pre-instability evolution of gi¬ 
ant planets. The map is displayed in the resonant (a r ,e r ) plane 
corresponding to h = 2.5°. Note that only the displacement in 
the resonant inclination 8I T is used to represent the dynamical 
stability. 

inwards, the resonance follows and the resonant semima¬ 
jor axis a T of bodies residing inside decreases. This in turn 
changes the resonant eccentricity e r because of the adiabatic 
invariance of resonant orbits. These changes are systematic 
and have different amplitudes for different orbits. Conse¬ 
quently, they should not be incorporated when calculating 
the distance d. On the other hand, our numerical runs sug¬ 
gest that the resonant inclination / r of test particles does 
not undergo substantial systematic changes under the influ¬ 
ence of migrating giant planets. If a large change in inclina¬ 
tion is registered, it usually means that a secular resonance 
affected the orbit. The resulting dynamical map thus rep¬ 
resents what we need - the location of regions crossed by 
secular resonances and the stable islands lying in between 
for the time period before Jupiter’s jump. 

Our choice of boundaries in the phase space to construct 
the dynamical map was a E (3.32, 3.45) AU and e E (0, 0.6). 
We created 40 x 40 grid in the (a, e) plane and covered it 
uniformly with six particles per each cell, assigning a fixed 
value of I = 2.5° to all of them. As already mentioned, 
the planetary migration is not able to strongly change the 
inclinations, that is why we map the region of low inclina¬ 
tions where the observed long-lived population dominates. 
The resulting pre-instability dynamical map is displayed in 
Fig. [10] 

We randomly distributed a group of 2, 000 test parti¬ 
cles over the identified stable islands (see Fig. ED with low 
inclinations I < 5° to set up a synthetic primordial popu¬ 
lation for our simulation. The angular osculating elements 
were chosen in such a way that condition © holds and the 
osculating elements are therefore identical to the resonant 
elements at t = 0. 

Evolution of test particles. 2Myr after the onset of the 
migration, the orbits are not significantly dispersed which 
is in agreement with our aim to study only the dynamically 
stable resonant bodies. After 6.1 Myr, the close encounters 
of giant planets begin to occur and start to perturb the 
islands inside the J2/1. Progression of Jupiter’s jump par¬ 


Ainit (D > 5 km) 

Armodel 

iV A 

7\rmodel 

iV B 

2000 

0 

0-1 

5000 

0 

1-2 

10000 

0 

2-4 

100000 

1 - 3 

18 - 35 


Table 2. The numbers N™° del a nd ]\f^ odel 0 f the asteroids sur¬ 
viving in the islands A and B up to the present as predicted 
by our model of the primordial population survival. The table 
represents how the resulting population changes with increas¬ 
ing number Aj n j ti of initial primordial asteroids. The bodies with 
D ^ 5 km are considered. First line roughly corresponds to the 
primordial population with particle density of the present outer 
main belt, the next two cases approximately con sider the particle 
density proposed in lMinton fc Malhotral lt2010lf . In the last line, 
we assume particle density of primordial main bel t shortly after 
its creation, as suggested bv lMorbidelli et~a~l . $2009); we note that 
this case is not very probable. 


tially destabilizes the islands and enables significant number 
of bodies to leak out and also to populate other regions of 
the 2:1 resonance. At the end of migration, 85 per cent of 
the initial population is lost and the rest is dispersed all over 
the resonant region, except for the low eccentricity region. 
This indicates that the only possible change in resonant ec¬ 
centricity during the migration of Jupiter is to increase. 

Surviving population. Following our procedure from Sec¬ 
tion E2 we used the same dynamical map to identify the 
long-lived asteroids (see Fig. ED and we applied the rescal¬ 
ing and the long-term dynamical decay. Because the parti¬ 
cle density in a hypothetical primordial population is not 
known, we varied the number of initial synthetic long-lived 
particles and calculated the expected number of asteroids in 
islands A and B surviving up to the present. 

The selected results are listed in Tableland the com¬ 
parison is made for asteroids with D > 5 km (we remind the 
reader that there are Nj^ s = 2 and Ng bs = 71 larger than 
5 km in the obs erved population). If we as sume the parti¬ 
cle density from iMinton fe Malhotral (l2010l ). which we also 
used in Section [372] the contribution of the primordial aster¬ 
oids to the observed population is negligible after long-term 
dynamical evolution. 

In order to observe a substantial contribution, the 
particle density in the J2/1 would have to be consid¬ 
erably larger than in the neighbouring main belt (at 
least ten times larger). Such particle density gradient 
is not very probable, because the 2:1 resonance would 
have to be dynamically protected against depleting mech¬ 
anisms and perturbations, arisi ng e.g. fro m planetary 
embryos (lO’Brien, Morbidelli fc Bottkef 120071 ') or due to 
the compact primordial configuration of planetar y orbits 
llMasset fe Snellgrovell200ll ; iRoig fe Nesvornvl|2014l ). 

We also investigated if A-island orbits can become B- 
island and vice versa. We found that < 1 per cent of the 
asteroids initially placed inside island A survive there and 

1 per cent drift into B island during the migration. On the 
other hand, the primordial population of island B preserves 

2 per cent of the original bodies and 2 per cent of them 
populate island A. Because the primordial B island is larger 
than the A island, it harbours more primordial asteroids 
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Figure 11 . A simulation of the primordial population survival in the fifth giant planet scenario llNesvornv fc Morbidellill2012h . The 
evolution of primordial long-lived orbits in the (a r ,e r ) plane is displayed. The vertical line indicates the migrating 2:1 resonance. The 
initial population was placed inside the islands identified in Fig. 1101 top left panel represents its orbital distribution 2 Myr after the onset 
of the migration. The following panels show the situation during Jupiter’s jump and also the final state. 


under the assumption of homogeneous particle density. The 
contribution of island B to the surviving long-lived popula¬ 
tion therefore dominates. Also note that after the migration 
the ratio Na/Nb — 1, thus the migration itself is not able to 
create an asymmetric population out of surviving primordial 
asteroids. Further orbital evolution is needed in this case. 

The orbital distribution in the (a r , 7 r ) plane at the end 
of the planetary migration indeed resembles the initial in¬ 
terval of inclinations, as both islands are populated by low -I 
orbits. To obtain high-/ orbits, which are observed in island 
A, we would have to assume the presence of their analogues 
in the pre-migration islands. But this problem is redundant 
since the primordial population probably do not survive up 
to the present, as we argued above. 


4 EFFECTS OF JUPITER-SATURN GREAT 
INEQUALITY 

In Section [3J we investigated evolution of the resonant pop¬ 
ulation during a major instability of the planetary system 
known as Jupiter’s jump. In the original mi grati on experi¬ 
ments performed in lNesvornv fc Morbidellil ( 2012ll . the vio¬ 
lent evolution of planetary orbits was usually followed by a 
residual smooth migration during which giant planets slowly 
approached their current orbits. We are in a similar situa¬ 
tion, the planetary configuration at the end of our insta¬ 
bility simulations slightly differs from the observed one. In 
principle, this fact does not affect our ability to locate and 


describe the long-lived resonant population of test particles 
because we employed fast and efficient method of dynamical 
mapping. The natural question then arises, whether the late 
migration stages can invoke perturbations that would affect 
the overall stability of the 2:1 resonance significantly. 

In particular, we want to check whether the period 
P a of libration of the resonant asteroids can be compa¬ 
rable with the Jupiter-Saturn great inequality period Pgi 
(which is the period of circulation of the 2Aj — 5As angle) 
and what effect would it have on the resonant population 
llFerraz-Mello et al.lll998l . see also Section [TJ). 

In order to mimic the smooth late planetary migration 
and lead giant planets towards their current orbits, we used 
a modified version of the SWIFT_rmvs3 integrator developed 
in iBroz et~akl J201lll . The integrator introduces an ad hoc 
dissipation term which modifies planetary velocity vectors 
v in each time step At according to the following relation: 


v (t + At) 


v(t) 


Av At 

1 -I-exp 

V Tmig 


Tmig ) 


(7) 


where Av = GM /ai n it — \JGM /afl n is the total dissipation 
determined as the difference of the initial and final mean ve¬ 
locity, T m ig is the migration time-scale, t is the time variable 
and to is an arbitrary initial time. The eccentricity damp¬ 
ing is also included and can be set independently for each 
planet by ch oosin g the damping parameter denoted as edamp 
rtMorbidelli et al.ll2010l ') . 

To set up the model, we used the final configuration 
of giant planets and test particles from our simulations of 
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Figure 12. Top: The same post-migration dynamical map as 
in Fig. [9] The primordial test particles surviving the jumping- 
Jupiter instability are plotted over the map. Circles represent 
candidates for long-lived orbits and crosses are the surrounding 
orbits. Bottom: the resonant semimajor axis a r vs sine of the 
resonant inclination sin I T . We plot only the particles identified 
above as long-lived. Black circles indicate those residing in the 
island B and gray open squares indicate those residing in the A 
island. 
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resonant capture^, but we only selected test particles lo¬ 
cated in broader surroundings of the stable islands to dis¬ 
card major part of short-lived asteroids. After this proce¬ 
dure, we launched a set of integrations with different val¬ 
ues of edamp and different migration time-scales, bearing 
in mind that the time-scale of the original experiments in 
iNesvornv fe Morbidellil (l2012| j was r m i g ~ 30 Myr. When the 
integrations finished, we selected only the runs in which gi¬ 
ant planets ended up with orbital parameters similar to the 
observed ones and we investigated the results at the time 
tfin when Pqi = 880 yr and the ratio of the Saturn’s and 
Jupiter’s orbital periods is approximately Ps/Pj — 2.49. 

In the following, we will discuss results of two runs with 
tfi n — 27 Myr and tfl n — 58 Myr. Fig. [13] shows the temporal 
evolution of the number N of test particles in these two runs, 
the GI period Pgi and the libration period P a of two typi- 


5 The results of simulations with primordial asteroids are not 
considered here, as these asteroids probably do not significantly 
contribute to the observed population. 



Figure 13. Results of two simulations of smooth late migra¬ 
tion and its effects on the resonant population. Two runs with 
tfin ~ 27 Myr (red and black curves and arrows) and t (\ n cs 58 Myr 
(blue and gray curves and arrows) are presented (migration time- 
scales T m ; g = 30 Myr and r m ; g = 60 Myr were initially chosen 
but the integration is evaluated only until the GI period reaches 
880yr). Top: temporal evolution of the number N of test parti¬ 
cles. Bottom: temporal evolution of the GI period Pqi and the 
libration period Pa. The latter is plotted for two representative 
cases. The arrows in the top panel approximately mark the time 
period during which Pqi — Pa- We note that ‘real’ Pqi evolves 
in an oscillatory manner but only with a small amplitude; a poly¬ 
nomial fit of the real Pqi evolution is plotted here for clarity. 


cal resonant asteroids. While Pa oscillates around a nearly 
constant value ~ 420 yr, Pqi initially lies below this value 
and rises smoothly as Jupiter and Saturn undergo divergent 
migration, and their eccentricities are damped. 

The evolution of the N(t) curves does not follow any 
simple decay law. This may serve as a proof that the level of 
chaotic diffusion inside the islands indeed evolves during the 
smooth late migration. The steepest slope of the N(t) depen¬ 
dency is reached during the time interval when Pgi — Pa- 
But one can also notice that the slope changes are smooth 
rather than sharp; the destabilization is not strictly bounded 
by the P gi — Pa condition. This is in agreement with the 
results of Ferraz-Mello et all (Il998l l who argued that even a 
near-resonant state of Pgi and P a can enhance the chaotic 
diffusion. However, it is clear that the number of asteroids 
surviving the smooth late migration depends mainly on the 
length of the time interval during which Pgi — P CT ; the pop¬ 
ulation is more depleted in the run with tfi n ~ 58 Myr. 

With the aim to evaluate the result of the models with 
residual migration included, we first identified test particles 
which were located inside the stable islands at the time tfi n . 
At this time, the orbital architecture of giant planets ap¬ 
proximately corresponds to the observed one and therefore 
the identification of long-lived asteroids can be achieved us¬ 
ing the current dynamical map. Finally, we derived expected 
numbers J\[^ odel an( l N™ odel 0 f asteroids that should be ob- 
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servable in the stable islands after 4 Gyr of orbital evolution. 
For this purpose, we rescaled the initial population of test 
particles and we applied simple exponential dynamical decay 
law, exactly as in Section [3.21 

We calculated following values for asteroids larger 
than 5 km: N£ oAel = 0-1 and Ng odel = 40-78 in case of 
t fln ~ 27 Myr and JV£ lodel = 0-1 and ATg lodel = 22-42 in case 
of tfln — 58 Myr (while the observed values are N^ s = 2 and 
77g bs = 71). We can conclude that the hypothesis of resonant 
capture is still valid and explains the existence^ of the long- 
lived population, if the giant planets finished their late stage 
migration on a time-scale comparable to r m ig — 30 Myr. 
Larger time-scales lead to a slower passage of the GI period 
over the values of the J2/1 libration period and this in turn 
causes an enhanced depletion of the stable islands. Even in 
such a case, a significant part of the observed population 
should originate from the resonant capture. 


5 COLLISION AL MODELS 

Our results of Section [3] and [I] imply that the long-lived 
J2/1 population was created by resonant capture during 
planetary migration. In this section, we further develop the 
framework of this hypothesis. We study collisional evolution 
of the J2/1 in order to confirm whether the size-frequency 
distribution of the long-lived resonant asteroids can survive 
a time period of 4 Gyr in a non-stationary state. We thus 
have to also accou nt for the epoch of the late heavy bom - 
bardment (LHB) dGomes et al.l 120051 : iLevison et al] l2009l f 
during which the transneptunian disc was destabilized and 
the flux of cometary projectiles through the solar system 
increased substantially. 

In the following collisional models, three populations 
of minor solar-system bodies are included: main-belt aster¬ 
oids, transneptunian comets and Zhongguos, i.e. we consider 
only a subset of the long-lived population. The reason for 
this simplification is that the observed SFDs of long-lived 
population and its dynamical subgroups share similar fea¬ 
tures (see Fig. [I]), but Zhongguos have the steepest slope 
7 = —5.1. Solving an inverse problem, we aim to explain 
the formation of an initial SFD with an even steeper slope, 
which would evolve towards the observed one due to colli¬ 
sions over the 4 Gyr timespan. Because of certain level of 
parametric freedom in collisional models, the same process 
of formation should also apply to Griquas with shallower 
observed SFD. 

5.1 Intrinsic probabilities and impact velocities 

As the first step in the construction of a collisional model, we 
compute the intrinsic probability Pi and the mean impact 
velocity Vi of colliding main-belt and long-lived reson ant as¬ 
teroids. We adopt the method of lBottke et al.l (1 199414 based 
on a geometrical form alism of orbital encounters introduced 
by iGreenberd (1 1982tl . To make our samples of orbits large 

6 Although the values derived here do not overlap the observa¬ 
tions in case of island A, the difference is only one asteroid. We 
think that this discrepancy is not significant as we are comparing 
small numbers and the evolution is definitely stochastic. 


colliding 

populations 

p 

[10—1-8 k m - 2 y r -!] 

Vi 

[kms -1 ] 

Zhongguos vs MB 

3.82 

5.36 

Griquas (I ^8°) vs MB 

3.65 

5.52 

Griquas (I >8°) vs MB 

1.81 

7.57 

average 

3.09 

6.15 


Table 3. The intrinsic probability P; and the mean impact ve¬ 
locity Vi for collisions of main-belt and long-lived J2/1 asteroids. 
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f/ Myr 

Figure 14. Temporal evolution of the intrinsic probability P;(f) 
(bottom curve) and the mean impact velocity Vj(t) (top curve) 
for collisions between transnept unian comets an d main-belt as¬ 
teroids, as it was calculated in lBroz et al.l JlO I :] ). We emphasize 
that we modify these dependences by numerical factors to allow 
for realistic lifetimes of comets. The values P;/3 and Vi/1.5 serve 
as an input for our models. 

enough, we input all long-lived orbits and the first 50, 000 
main-belt objects from the AstOrb catalogue. 

We split the long-lived population into Zhongguos, Gri¬ 
quas with inclination 1^8° and Griquas with inclina¬ 
tion I > 8° and calculate the intrinsic probability Pi and 
weighted mean impact velocity Vi for them, separately. Us¬ 
ing this separation, we want to check if the intrinsic proba¬ 
bilities for Griquas may differ from Zhongguos significantly. 

The resu lts are s ummarized in the Table [3] For refer¬ 
ence, [Dahlgren| (119981') computed Pi = 3.1 x 10 -18 km -2 yr -1 
and Vi = 5.28kms _1 for collisions between main-belt as¬ 
teroids. We can see that both Zhongguos and low-inclined 
Griquas have Pi and Vi only slightly higher than these refer¬ 
ence values. This is caused by the moderate values of eccen¬ 
tricities in the J2/1 population. The orbits then more likely 
intersect with those in the main belt. The higher collisional 
velocity is also plausible because we are comparing an outer 
main-belt population with the rest of the main belt. 

On the other hand, Griquas with high inclinations can 
avoid intersecting some of main-belt orbits, thus their in¬ 
trinsic probability is lower by about a factor of two. Even 
at this point, we can conclude that the observed difference 
in slopes of SFDs (when Griquas and Zhongguos are sepa¬ 
rated) cannot be explained by the differences in Pi and V 
because the shallower SFD of Griquas would require higher 
Pi in order to collide more often. 
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Figure 15. Results of collisional models created with the Boul¬ 
der code, plotted as temporal evolution of the SFDs. Description 
of individual curves is given in the legend of each plot. The evolv¬ 
ing curves are displayed in the time of the simulation t , which is 
also shown. Top: synthetic Zhongguos are created by capture from 
the main-belt background population affected by collisions with 
transneptunian comets (see also Section I5.2D . Middle: the situa¬ 
tion from top left panel is reproduced, but a single D ~ 100 km 
asteroid is added to Zhongguos and we study the possibility of its 
catastrophic breakup over 4Gyr (see also Section |5.3|> . Bottom: 
we demonstrate what parameters of the initial SFD of Zhongguos 
are needed in order to evolve it towards the observed one due to 
collisions (see Section l5.4i for discussion). 
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5.2 Capture from the main-belt population 

Let us ask a question whether the J2/1 SFD can originate 
from the main-belt SFD. We assume that the capture of 
the resonant population is not size dependent and thus the 
captured SFD resembles that of the main belt at the time 
of planetary migration, scaled down by a numerical factor. 
Our arbitrary choice of this factor is such that the number 
of the largest captured bodies is approximately the same as 
we observe in the J2/1 population. This allows us to imme¬ 
diately compare the slope of captured and observed SFD of 
Zhongguos. 

The main belt itself was likely more populous 4Gyr 
ago (b y a factor of three according to iMinton fc Malhotral 
d2010f) 1 but has r emained in a near-s tationary collisional 
regime ever since (Bott ke et al.ll200 5ah The only possible 
stage of evolution, which could have temporarily increased 
the slope of the main belt SFD, was the late heavy bombard¬ 
ment. We therefore investigate this early period of collisions 
between the main belt and transneptunian comets. We fo¬ 
cus on diameters D < 25 km to see whether the steepness 
in this interval can be instantaneously increased - if so, a 
population captured from such SF D would adopt that slo pe. 

We use the Boulder code dMorbidelli et al] 120091 ') to 
construct a suitable collisional model. We setup initial SFDs 
by a piecewise power law function which is characterized by 
three differential slope indices g a , qb and q c in the intervals 
of diameters D > D\, (D 2 ,Di) and D < D 2 , respectively. 
The SFD is normalized by setting the number N noiIn of bod¬ 
ies larger than D\. The following values are used in case of 
the main-belt SFD: D\ = 100 km, D 2 = 14 km, g a = —5.0, 
qb = —2.3, q c = —3.5, Nnorm = 1110. A few D ~ 1000 km as¬ 
teroids are added in order to properly reprod uce the present 
state. Using the same assumptions as lBroz et al.l (l2013ll . the 
cometary disc is characterized as: D 1 = 100 km, q A = —5.0, 
qb = q c = -3.0, N n orm = 5 x 10 7 . 

Because the evolution of the transneptunian disc is 
dominated by its fast dynamical dispersion, the intrinsic 
probability -Pi(i) and mean impact velocity Ui(f) of collisions 
with comets are time-dependent quantities. Th e tempora l 
evolut ion of -Pi(i) and V5(f), which was d e rived in lBroz et al.l 
(l2013h . is shown in Fig. 1141 iBroz et~akl (l2013l l also argued 
that it is necessary to modify this dependence in order to 
mimic the effects of spontaneous cometary breakups due to 
various physical processes. Regarding these effects, they de¬ 
rived quantities P(t) = Pi(t)/3 and Vi(t) = Vi(i)/1.5 as a 
feasible modification. 

The specific energy Q£, required to disperse 50 per cent 
of th e shattered target is des cribed by the polynomial scaling 
law (iBenz fe Asphauglll999l l 

Qd (r) = — (Qor“ + Bpr b ) , (8) 

(/fact ' ' 

where r denotes the target’s radius in cm. Material param- 
eters p, Q 0 , a, B, b and Qf ac t for basaltic rock (used for 
asteroids) and water ice (used for comets) are listed in Ta¬ 
bic [4] 

We simulated the collisional evolution for a time period 
of 10 Myr in which Pi (t) reaches its maximum (see Fig. [14] 
and top panel of Fig. m . The flux of cometary projectiles in¬ 
duced by ongoing planetary migration reaches its climax and 
strongly affects the main-belt SFD. We performed 100 nu¬ 
merical realisations of the model with different seeds of the 
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p 

[gem -3 ] 

Qo 

[ergg- 1 ] 

a 

B 

[ergg -1 ] 

b 

<7fact 

D pb 

[km] 

^disrupt 

[km] 

Aproject 

N co i 

Basalt: 






50 

6 

25821 

0.18 

3.0 

7 x 10 7 

-0.45 

2.1 

1.19 

1.0 

70 

9 

7837 

0.10 

Water ice: 






100 

15 

2589 

0.07 

1.0 

1.6 x 10 7 

-0.39 

1.2 

1.26 

3.0 

120 

19 

1777 

0.07 

Table 4. Material parameters p, 

Qo, B, b and <j fact of 

Table 5 

. The results of a stationary collisional model: the number 

the polynomial 

scaling 

law (see equation JHJ) adopted from 

N co i of catastrophic breakups of a parent body with the diameter 

Benz &; Asohaug d 19991). 

The first line is used in case of main- 

Dpb over 4 Gyr timespan. We also list the minimal size disrupt 


belt or resonant asteroids and the second line is used for comets. 


and the number .V pro j ecIi of suitable projectiles. 


built-in random-number generator to account for stochastic- 
ity of the evolution and possible low-probability breakups. 
At the end of the simulation, we assume the resonant cap¬ 
ture takes place with the efficiency factor 10 -4 . The steep 
part of the captured SFD can be characterized by the slope 
7 = —3.0. Although the region of the km-sized bodies is 
steeper with respect to the initial state of the main-belt 
SFD, it is not steep enough to reach the slope of the SFD of 
observed Zhongguos. We thus consider this scenario unlikely. 

5.3 Catastrophic disruption of a captured parent 
body 

From the result of the previous section, it is obvious that 
a different explanation of the steep initial SFD is needed. 
Here we test the possibility that a large parent body was 
captured inside the 2:1 resonance by chance, it was subse¬ 
quently disrupted and its fragments formed a family. 

We first use a stationary model to investigate the prob¬ 
ability of such a catastrophic breakup event. As a prereq¬ 
uisite, we estimated a lower limit of the parent body size 
Dpb. To achieve that, we searched th e outcomes o f hy- 
drodynamic disruption models in iDurda et al.l ll2007l ') and 
iBenavidez et al.1 ll2012l f for SFDs which have approximately 
the same slope as the observed SFD of Zhongguos. We 
rescaled the sizes of synthetic asteroids in these datasets so 
that the diameter of the largest remnant would correspond 
approximately to the size of (3789) Zhongguo. We selected 
all reasonable fits and derived range of admissible values 
Dpb £ (50,120) km with the median value Dpb = 70 km. 

We then use the following relation expressing the num¬ 
ber A co i of catastrophic disruptions of parent bodies with 
the diameter Dpb due to collisions with a population of 
projectiles: 

D 2 

AT col = PilVpBiVproject -p At, (9) 

where P denotes the intrinsic probability, Npb is the num¬ 
ber of parent bodies, Ap ro j e ct is the number of projectiles 
capable of disrupting the parent body and At is the consid¬ 
ered timespan. 

Since there are no large bodies observed in the cur¬ 
rent J2/1 population, we assume IVpb = 1. To provide the 
lower limit o n the diameter of the projectiles (Disrupt, we use 
the model of lBottke et al.l d2005bh in com bi nation with th e 
scaling law given by equation <( 8 j) dBenz fe Asphau dll999lh 
We then take -/Vp ro j ec t for different values of (Disrupt as the 
number of main-belt asteroids with diameters D ^ (Disrupt- 
Several resulting values of N co \ over At = 4Gyr timespan 


are given in Tabic [5] It turns out that the probability of a 
catastrophic disruption is ^ 18 per cent when considering a 
single target only. 

Note that a more realistic case should take into consid¬ 
eration that there is no observable evidence of a collisional 
cluster, thus the hypothetical break u p event m ust have oc¬ 
curred more than 1 Gyr ago (IBroz et aid 1200511 . Of course, 
one should also account for the influence of cometary flux 
during LHB which can temporarily increase the rate of col¬ 
lisions but also tends to speed up the evolution of SFDs. 

We therefore test the possibility of a catastrophic 
breakup once again in the more sophisticated framework of 
the Boulder code and we also check the influence of comets. 
The setup for main-belt asteroids and comets is the same as 
in Section E21 as well as the SFD of synthetic Zhongguos, 
which is only modified by adding a single D ~ 100 km aster¬ 
oid. The intrinsic probability and the mean impact velocity 
for MB vs. Zhongguos collisions are taken as the average 
of values from Table [3] We simulated 4 Gyr of collisional 
evolution. 

The middle panel of Fig. 1151 shows a set of 100 realisa¬ 
tions of our model. At the beginning of the simulations, the 
synthetic population of Zhongguos is strongly affected by the 
cometary flux, which mostly causes cratering of the large 
body in the resonance. Moreover, a family-forming event 
takes place in a few runs. Simultaneously, the comets speed 
up the collisional evolution in and below the region of mid 
sized asteroids. As a result, families inside the J2/1 are com- 
minutioned too fast and the steepness of their SFDs drops 
below the observed one. 

When studying the family-forming events after the 
LHB, we focused only on cases with the size of the largest 
fragment (or remnant) Dlf SS 10 km to obtain a body simi¬ 
lar to (3789) Zhongguo. In three such cases, which occurred 
at the time of simulation t = 2.1, 2.8 and 3.1 Gyr, a parent 
body with the diameter Dpb — 85 km was shattered into a 
family with the largest fragment having Dlf — 18 km and 
the largest remnant having Dlr — 60 km. Corresponding 
evolved SFDs matclQ the observed one, except for the pres¬ 
ence of the parent-body remnant. 

We can conclude that a family-forming event inside the 
2:1 resonance is a very unlikely process because a complex 
set of constraints has to be satisfied: A major collision would 
have to occur despite its low probability, which is typically 

7 In fact, the SFDs resulting from our simulations lie slightly be¬ 
low the observed case but this discrepancy can be easily removed 
assuming a bit larger parent body. 
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Figure 16. Dependence of the absolute magnitude H on the 
proper semi-major axis a p for a synthetic family (top panel). The 
family is used in a simple Monte-Carlo test, in which we study 
how the shape of the SFD changes if we select only a part of the 
family. An example, for the part of the family bordered by the 
vertical lines, is given in the bottom panel. We compare cumu¬ 
lative size-frequency distributions of the entire synthetic family 
(dashed curve) and of its part (solid curve). The steep segments 
are approximated with a power law function and obtained slope 
indices are shown for reference. 

10 per cent, as given by the stationary model, or 3 per cent, 
as derived from the more sophisticated model. The created 
SFD would have to be very steep, preferably without pres¬ 
ence of a parent-body remnant. If the remnant was present, 
one would have to rely on its subsequent elimination due to 
dynamical depletion because there is no large asteroid ob¬ 
served in the J2/1. Finally, the breakup would have to take 
place after the LHB, but early enough to allow for dispersion 
of the collisional cluster. 

5.4 Capture from a family 

Previous collisional models imply that it is difficult to ex¬ 
plain the steep SFD of Zhongguos by standard processes 
(such as direct capture from the background main-belt pop¬ 
ulation or family-forming event) that could have occurred 
during resonant capture or later on. We thus employ a sort 
of ‘reversed’ method in this section. We first try to discover 
what parameters of the initial SFD are needed in order to 
obtain the observed one after 4 Gyr of collisional evolution. 


Then we discuss possible process which could have led to 
formation of such an SFD. 

Using trial and error, we established the following ini¬ 
tial SFD as an appropriate one: D\ = 10 km, D 2 = 5 km, 
g a = —6.6, <?b = —6.6, q c = —3.5, = 12. The change of 

slope at D 2 is necessary to reasonably satisfy the mass con¬ 
servation law. The main-belt and cometary SFDs are the 
same as in the previous models. 

The result of 100 runs of the collisional code is displayed 
in bottom panel of Fig. [To] The evolved SFD of the resonant 
population corresponds to the observed one very well. The 
steep parts have the same slope, the only difference can be 
seen in the region of small bodies where the synthetic SFD 
is abundant. But one has to realize that the observed SFD 
is probably biased in this region. 

The mechanism responsible for creation of the high ini¬ 
tial slope m ay be th e foll owing (and it is r elated to the effect 
described in lCarruba, Alibaae fe Souamil (120141) ). Let us as¬ 
sume resonant capture from a hypothetical outer main-belt 
family instead of capture from the background main-belt 
population. The stable islands are able to capture only a 
fraction of this family because jump of the resonance is 
thought to be fast (it does not sweep through the resonance), 
the capture efficiency is limited (see Section [3.21) . and also 
the width of the islands in semimajor axis is relatively small 
(see e.g. Fig. [T|). If the captured part of the family is located 
farther away from the original position of the parent body, 
its SFD will be very steep because smaller fragments have 
higher ejection velocities and therefore land at orbits that 
are more distant from the parent body. 

We employed a simple Monte-Carlo test to investigate 
the described possibility. We generated a uniform distribu¬ 
tion of 10, 000 test particles satisfying the bounds 2.52 ^ 
a ^ 2.70 AU and log (|a — a c | /2 x 10 -4 ) /0.2 ^ H ^ 16 in 
the (a, H) plane where a c — 2.64 AU is value of the central 
semimajor axis and H is the absolute magnitude. Finally, 
we assigned diameters D to the test particles on the basis 
of the H distribution, assuming single value of the geomet¬ 
ric albedo pv = 0.05 throughout our sample. This way we 
created a synthetic family (similar e.g. to Eunomia family). 
Then we randomly moved a A a = 0.02 AU window in the 
semimajor axis over the collisional cluster and we monitored 
the SFDs in the successively selected regions. 

An example is given in Fig. 1161 Obviously, extremely 
high slope can be reached by this process. The resulting 
slope is well above the estimated lower limit which is needed 
for the initial SFD of Zhongguos. Note that this mechanism 
does not require the family itself to have steep SFD. The 
steepness is achieved afterwards by selective resonant cap¬ 
ture from an appropriate region. 

Finally, we remark that an interesting link can be found 
between the hypothesis of capture from a family and our dy¬ 
namical simulations. In Section we argued that a source 
population for resonant capture should contain greater num¬ 
ber of asteroids with low inclination in order to explain the 
observed high concentration of low -I B-island bodies. Thus 
we can conclude that if the hypothetical outer main-belt 
family was indeed captured, it probably must have been lo¬ 
cated on low inclinations. 
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6 CONCLUSIONS 

Let us briefly review the content of this paper. We updated 
the population of asteroids residing in the 2:1 mean-motion 
resonance with Jupiter using recent observational data. The 
new list of resonant objects now contains 370 bodies, which 
can be divided on the basis of their mean dynamical lifetime 
to 140 short-lived and 230 long-lived asteroids. Our revision 
of physical properties of the resonant population is generally 
in agreement with the conclusions of previous studies; new 
result is our estimate of the mean albedo py = (0.08 ±0.03), 
based on the data from the WISE database. 

The long-term dynamics of the two quasi-stable islands 
A and B was studied, using the results of Skoulidou et al. (in 
preparation) who took into account also the perturbations 
induced by the terrestrial planets and the semimajor axis 
drift caused by the Yarkovsky effect. The population in the 
islands decays exponentially, but the escape rate is signifi¬ 
cantly faster in island A, the e-folding times being 0.57 Gyr 
for island A and 0.94 Gyr for island B. Hence, the dynami¬ 
cal evolution on Gyr-long time-scales results into differential 
depletion of the islands, which is certainly one of the reasons 
for the observed asymmetric A/B population ratio. 

The primary goal of this paper was to explain the for¬ 
mation of the long-lived population, satisfying the observa¬ 
tional constraints. We tested two hypotheses: 

(i) capture from the outer main belt, 

(ii) survival of primordial long-lived resonant asteroids, 

both i n the framewo rk of the fi ve giant planets migration sce¬ 
nario (INesvornv fc Morbidellil[2012l ) with a jumping-Jupiter 
instability. 

Our simulations of the instability imply that both pro¬ 
cesses could have been at work. The capture itself, to¬ 
gether with subsequent differential depletion due to long¬ 
term dynamical evolution, can explain the observed pop¬ 
ulation in both islands, if we assume that the number of 
asteroids in the main belt 4 Gyr ago w as three times larger 
than the current one, as suggested by [Minton fe Malhotral 
(l20ldl . Namely, the numbers of asteroids with D ^ 5 km, 
which should survive to this day, are 7V™ odel = 1-3 and 
^ymodei _ g2_i21, as predicted by our model. For compari¬ 
son, the observed numbers are N^ s = 2 and IVg ba = 71. 

We also modeled the late residual planetary migration 
to check on the effect of possible resonance between the li- 
bration period in the 2:1 commensurability and the raising 
period of Jupiter-Saturn great inequality. The impact of this 
effect on the long-lived population depends on the time-scale 
of the late migration. If the time-sca le was ~ 30 Myr (which 
corres ponds with the scenarios in iNesvornv fe Morbidellil 
l)2012l) 1. the destabilization of the islands would be weak 
enough for the captured population to survive in a state 
similar to the observations. More specifically, if we account 
for the effect of great inequality in our model, the resulting 
values are JV“ odel = 0-1 and N™ odel = 40-78. If the time- 
scale was ~ 60Myr, the chaos in the stable islands would 
be intensified for longer period of time and the depletion of 
the captured population would be more significant, leading 
to N™ odel = 0-1 and Ag lodel = 22-42. 

In case of primordial resonant objects, we discovered 
that approximately 1 per cent of A-island and 4 per cent 
of B-island asteroids survive Jupiter’s jump on long-lived 


orbits. However, we argued that the primordial asteroids 
probably do not contribute to the observed population at 
all. The reason is that the primordial population would have 
to exhibit particle d e nsity larger than the one proposed by 
iMinton fe Malhotral (hold ) by a factor of ten, in order to 
survive 4 Gyr of post-migration dynamical evolution. Such 
particle density would create a gradient with respect to the 
neighbouring main belt which we assume dubious. 

Finally, by creating several collisional models, we 
demonstrated that the observed steep SFD of the long-lived 
asteroids cannot be explained by capture from the back¬ 
ground main-belt SFD, affected by the late heavy bombard¬ 
ment. We also proved that a family-forming catastrophic 
disruption inside the J2/1 is very unlikely. 

Our main conclusion is that the long-lived J2/1 popu¬ 
lation was probably formed by capture from a hypothetical 
outer main-belt family during Jupiter’s jump. If this is the 
case, then the long-lived asteroids in the 2:1 resonance with 
Jupiter represent the oldest identifiable remnants of a main- 
belt asteroidal family. 

There are several improvements that are needed to con¬ 
clude the debate on the long-lived population in the 2:1 res¬ 
onance and its origin. For example, it is appropriate to as¬ 
sess the possibility of the Themis family formation event as 
a contributor to the resonant population. Although we did 
not study this hypothesis in detail, we summarized several 
results of our preliminary tests in Appendix [A] Our simu¬ 
lations suggest that this possibility is not viable and this in 
turn supports the hypothesis of resonant capture. 

Other improvements of our work might include a devel¬ 
opment of a self-consistent model of both dynamical and col¬ 
lisional evolution of long-lived asteroids. Finally, we suggest 
to study the hypothetical primordial population in the epoch 
before the jumping-Jupiter instability in order to properly 
estimate primordial particle density inside the stable islands. 
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APPENDIX A: CONTRIBUTION OF THE 
THEMIS FAMILY EJECTA TO THE 
LONG-LIVED POPULATION 

We shall briefly discuss a possible role of neighbouring 
Themis family in creation of the l ong-lived resonant pop¬ 
ulation. Although iBroz et al] (1200514 demonstrated that the 
fragments from Themis family cannot be transported to the 
stable islands by the Yarkovsky semimajor axis drift, a possi¬ 
bility of direct injection during the Themis family formation 
event is still an open question. During the formation, frag¬ 
ments ejected at high velocities (u e j — 400ms - , or more) 
may fall directly in the J2/1 resonance and thus contribute 
somehow to the stable population. There are several caveats, 
however. 

(i) One has to assume that a substantial part of 
fragments (~ 10 per cent) have very large ejection ve¬ 
locities with respect to (24) Themis, or the respective 
parent body with Dpb — 400 km. We tried to use 
an ‘extreme’ size-i ndependent velocity field prescribed by 
iFarinella, Gonczi fe Froeschle! d 19931) relation, with param¬ 
eters Vesc = 170ms -1 and a = 3.25. Then up to N(D > 
5 km) = 10-30 fragments land within island B, according 
to our tests. This number should be further decreased by a 
subsequent long-ter m orbi tal evolution, as Themis family is 
(2.5 ± 1.0) Gyr old llBroz et al.ll2013h . Let us also note that 
the collisional evolution of resonant bodies was accounted for 
automatically, as we did this simulation with the currently 
observed size-frequency distribution of Themis family. 

(ii) At the same time, it is required that the breakup 
takes place when the true anomaly /i mp ~ 0° for this ejection 
scenario to work; otherwise, the number of objects landing 
in the islands decreases as well. But this particular impact 
geometry does not seem to be compatible with the observed 
shape of Themis family, namely with a large eccentricity 
dispersion of the family below the J11/5 resonance, at a p = 
3.03 AU (see Fig.lATl). 

(iii) Starting with a more reasonable siz e -depen dent ve¬ 
locity field - as used in iBroz fe Morbidellil (I2OI3I) for Eos 
family, which has a similar parent body size - makes the 
contribution to the long-lived population as low as N(D > 
5 km) = 2. 
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resonant semimajor axis a r [AU] 

Figure Al. A simulation of a synthetic Themis-like family for¬ 
mation event. The figure represents the state 100 Myr after the 
breakup. The orbital distribution of observed Themis family is 
displayed in the proper elements (black squares), while the orbital 
distribution of the synthetic family is displayed in the resonant el¬ 
ements (blue circles). This causes the mutual shift in eccentricity. 
Note the differences in extent of both families and the eccentricity 
dispersion beyond the 11:5 mean-motion resonance with Jupiter. 


(iv) There is no way to explain the existence of aster¬ 
oids in island A by the ejection, both the eccentricities and 
inclinations of the observed A-island objects are too large. 
Therefore, the ejection from Themis is not a viable hypoth¬ 
esis in case of island A. 

This paper has been typeset from a T)eX/ I^T^X file prepared 
by the author. 








